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LONG-TERM  GOALS 

I  am  part  of  a  multi-institutional  research  team  funded  by  the  ONR-sponsored  Radiance  in  a  Dynamic 
Ocean  (RaDyO)  program.  The  primary  research  goals  of  the  program  are  to  (1)  examine  time- 
dependent  oceanic  radiance  distribution  in  relation  to  dynamic  surface  boundary  layer  (SBL) 
processes;  (2)  construct  a  radiance-based  SBL  model;  (3)  validate  the  model  with  field  observations; 
and  (4)  investigate  the  feasibility  of  inverting  the  model  to  yield  SBL  conditions.  Our  goals  are  to 
contribute  innovative  measurements,  analyses  and  models  of  the  sea  surface  roughness  at  length  scales 
as  small  as  a  millimeter.  This  characterization  includes  microscale  and  whitecap  breaking  waves. 

The  members  of  the  research  team  are 

Michael  Banner,  School  of  Mathematics,  UNSW,  Sydney,  Australia 

Johannes  Gemmrich,  Physics  and  Astronomy,  UVic,  Victoria,  Canada 

Russel  Morison,  School  of  Mathematics,  UNSW,  Sydney,  Australia 

Howard  Schultz,  Computer  Vision  Laboratory,  Computer  Science  Dept,  U.  Mass.,  Mass 

Christopher  Zappa,  Lamont  Doherty  Earth  Observatory,  Palisades,  NY 

OBJECTIVES 

Nonlinear  interfacial  roughness  elements  -  sharp  crested  waves,  breaking  waves  as  well  as  the  foam, 
subsurface  bubbles  and  spray  they  produce,  contribute  substantially  to  the  distortion  of  the  optical 
transmission  through  the  air-sea  interface.  These  common  surface  roughness  features  occur  on  a  wide 
range  of  length  scales,  from  the  dominant  sea  state  down  to  capillary  waves.  Wave  breaking  signatures 
range  from  large  whitecaps  with  their  residual  passive  foam,  down  to  the  ubiquitous  centimeter  scale 
microscale  breakers  that  do  not  entrain  air.  There  is  substantial  complexity  in  the  local  wind-driven  sea 
surface  roughness  microstructure,  including  very  steep  nonlinear  wavelets  and  breakers.  Traditional 
descriptors  of  sea  surface  roughness  are  scale-integrated  statistical  properties,  such  as  significant  wave 
height,  mean  squared  slope  (e.g.  Cox  and  Munk,  1954)  and  breaking  probability  (e.g.  Holthuijsen  and 
Herbers,  1986).  Subsequently,  spectral  characterisations  of  wave  height,  slope  and  curvature  have  been 
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measured,  providing  a  scale  resolution  into  Fourier  modes  for  these  geometrical  sea  roughness 
parameters.  More  recently,  measurements  of  whitecap  crest  length  spectral  density  (e.g.  Phillips  et  al, 
2001,  Gemmrich  et  al.,  2008)  and  microscale  breaker  crest  length  spectral  density  (e.g.  Jessup  and 
Phadnis,  2005)  have  been  reported. 

Our  effort  seeks  to  provide  a  more  comprehensive  description  of  the  physical  and  optical  roughness  of 
the  sea  surface.  We  achieve  this  through  the  analysis  of  our  suite  of  comprehensive  sea  surface 
roughness  observational  measurements  within  the  RADYO  field  program.  These  measurements  cover 
the  fundamental  optical  distortion  processes  associated  with  the  air-sea  interface.  In  our  data  analysis, 
and  complementary  collaborative  effort  with  RaDyO  modelers,  we  are  investigating  both  spectral  and 
phase-resolved  perspectives.  These  will  allow  refining  the  representation  of  surface  wave  distortion  in 
present  air-sea  interfacial  optical  transmission  models. 

APPROACH 

I  am  working  within  the  larger  team  (listed  above)  measuring  and  characterizing  the  surface  roughness. 
We  build  substantially  on  our  accumulated  expertise  in  sea  surface  processes  and  air-sea  interaction. 
This  team  contributed  the  following  components  to  the  primary  sea  surface  roughness  data  gathering 
effort  in  RaDyO: 

•  polarization  camera  measurements  of  the  sea  surface  slope  topography,  down  to  capillary  wave 
scales,  of  an  approximately  lm  x  lm  patch  of  the  sea  surface  (see  Figure  1),  captured  at  video 
rates.  [Schultz,  Zappa] 

•  co-located  and  synchronous  orthogonal  75  Hz  linear  scanning  laser  altimeter  data  to  provide 
spatio-temporal  properties  of  the  wave  height  field  (resolved  to  0(0. 5m)  wavelengths)  [Banner, 
Morison] 

•  high  resolution  video  imagery  to  record  whitecap  data  from  two  cameras,  close  range  and  broad 
field  [Gemmrich] 

•  fast  response,  infrared  imagery  to  quantify  properties  of  the  microscale  breakers,  and  surface  layer 
kinematics  and  vorticity  [Zappa] 

•  air-sea  flux  package  including  sonic  anemometer  to  characterize  the  near-surface  wind  speed  and 
wind  stress  [Zappa] 

The  team’s  data  analysis  effort  includes:  detailed  analyses  of  the  slope  field  topography,  including 
mean  square  slope,  skewness  and  kurtosis;  laser  altimeter  wave  height  and  large  scale  wave  slope  data; 
statistics  of  whitecap  properties,  as  well  as  statistical  distributions  of  whitecap  crest  length  density  in 
different  scale  bands  of  propagation  speed  and  similarly  for  the  microscale  breakers,  as  functions  of  the 
wind  speed/stress  and  the  underlying  dominant  sea  state.  Our  contributions  to  the  modeling  effort 
focus  on  using  RaDyO  data  to  refine  the  sea  surface  roughness  transfer  function.  This  includes  the 
representation  of  nonlinearity  and  breaking  surface  wave  effects  including  bubbles,  passive  foam, 
active  whitecap  cover  and  spray,  as  well  as  micro-breakers. 

WORK  COMPLETED 

My  effort  in  FY10  has  been  primarily  in  the  analysis  of  the  surface  video  data  obtained  during  the 
RaDyO  field  experiment  off  Hawaii,  August  28  -  September  16,  2009  as  well  as  refinement  to  the  data 
analysis  for  the  field  experiment  in  the  Santa  Barbara  channel  during  September  4-  28,  2008. 1  also 
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participated  in  the  ONR  Optics  review  symposium  (February  19-20,  2010)  and  the  RaDyO  workshop 
(June  1  -  3,  2010)  where  I  presented  a  preliminary  comparison  of  the  results  from  our  analysis  of  the 
RaDyO  -  Santa  Barbara  Channel  video  data  and  the  RaDyO  -  Hawaii  video  data. 

I  continued  work  on  the  analysis  techniques  for  characterizing  breaking  crests.  Routines  for  image 
rectification  and  the  extraction  of  length  and  propagation  speed  of  breaking  crests  have  been  adopted 
for  the  Hawaii  data  set.  Two  approaches  to  establishing  the  breaking  crest  length  distribution  were 
evaluated:  the  spectral  method  (Thomson  &  Jessup,  2009)  and  the  individual  tracking  method 
(Gemmrich  et  al,  2008).  All  final  processing  is  based  on  the  individual  tracking  method. 

Of  major  significance  to  our  group's  effort  was  the  deployment  of  our  polarimeter  in  the  RaDyO 
observational  periods  from  Scripps  Pier,  from  FLIP  in  the  Santa  Barbara  channel  and  off  Hawaii. 
Details  on  progress  with  this  development  are  given  in  the  companion  ONR  RaDyO  Annual  Reports 
by  Schultz  and  Zappa. 

RESULTS 

Figure  1  below  shows  a  schematic  of  the  instrumentation  deployed  in  the  field  experiments. 
Instrumentation  and  set-up  were  very  similar  in  the  Santa  Barbara  channel  (2008)  and  the  Hawaii 
(2009)  field  experiment.  Banner/Morison  deployed  two  orthogonal  line  scanning  lidars.  The  lidars 
were  positioned  on  the  boom  so  that  their  intersection  point  was  within  the  common  footprint  of  the 
polarimetric  (Schultz),  infrared  (Zappa)  and  visible  (Gemmrich)  imagery  cameras  which  were 
measuring  small-scale  surface  roughness  features  and  breaking  waves. 
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Figure  1.  Schematic  of  the  instrumentation  set-up  deployed  from  the  R/P  FLIP  starboard 
boom  during  the  Santa  Barbara  Channel  and  Hawaii  experiments.  The  end  of  the  boom  was 
about  9m  above  the  mean  water  level.  The  approximate  field  of  view  of  the  various  instruments 
is  shown.  A  second  wide  angle  whitecap  video  camera  was  mounted  on  the  crow’s  nest 
of  R/P  FLIP  approximately  26m  above  the  water  level  to  image  the  larger  whitecaps. 
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Zappa  deployed  his  infrared/visible  camera  system  and  his  environmental  monitoring  system  (sonic 
anemometer,  water  vapor  sensor,  relative  humidity/temperature  probe,  motion  package,  pyranometer 
and  pyrgeometer).  Gemmrich  deployed  2  video  visible  imagery  cameras.  One  camera  was  mounted  on 
the  main  boom  next  to  our  other  instrumentation  packages,  the  second  camera  was  mounted  higher  up 
to  view  larger  scale  breaking  events.  Schultz  deployed  an  instrument  package  located  on  the  boom  that 
includes  a  polarimetric  camera  imaging  the  very  small-scale  waves.  The  individual  data  acquisition 
systems  were  synchronized  and  the  various  data  sets  can  be  interrelated  to  within  0. 1  seconds.  In  this 
way  breaking  wave  properties  can  be  related  to  the  phase  of  the  underlying  waves. 

A  wide  range  of  conditions  prevailed  during  the  field  experiment  in  Santa  Barbara  channel  (September 
4-28,  2008)  where  the  wind  speed  Ufo  ranged  from  light  and  variable,  up  to  12  m/s.  The  scale  of  wave 
breaking  ranged  from  micro  breakers  to  small-scale  breakers  with  air  entrainment  to  breaking 
dominant  waves  with  up  to  2m  wave  height.  Environmental  conditions  during  the  Hawaii  experiment 
(Aug  28  -  Sep  16,  2009)  were  less  variable  with  a  nearly  fixed  wind  direction  and  wind  speeds  slowly 
fluctuating  between  8  m/s  and  12  m/s.  Periods  of  comparable  wind  speeds  resulted  in  significantly 
different  wave  fields  (significant  wave  height,  directionality,  dominant  period)  at  the  two  experimental 
sites. 

These  data  are  being  analysed  in  terms  of  breaking  crest  length  density,  foam  coverage  and  whitecap 
persistence. 

Preliminary  results  on  the  breaking  crest  length  distribution  A(c),  obtained  by  the  individual  tracking 
method,  are  shown  in  Figure  2.  All  curves,  for  both  data  sets,  show  the  maximum  distribution  to 
breaking  crests  in  the  intermediate  to  short  wave  range  (small  phase  speeds).  However,  the  absolute 
values  and  the  slope  of  the  curves  vary  significantly  throughout  the  experiments.  Generally,  the  slopes 
are  less  steep  than  the  canonical  value  A(c)  oc  A6  proposed  by  Phillips  (1985)  for  the  equilibrium 
range,  which  is  based  on  assuming  comparable  importance  of  wind  input,  nonlinear  spectral  transfer 
and  dissipation  through  breaking.  The  discrepancy  seems  to  be  somewhat  larger  for  the  coastal  wave 
field  in  Santa  Barbara  channel.  The  slope  decreases  with  increasing  wave  age,  i.e.  the  relative 
importance  of  large  scale  breakers  increases  as  the  wave  field  develops.  Thomson  et  al  (2009)  found 
the  A(c)-distribution  to  stratify  according  to  the  dominant  steepness  of  the  wave  field  kpHs .  Here,  we 

do  not  find  such  a  dependency. 

These  breaking  crest  length  distributions  are  solely  based  on  video  data.  They  are  also  compared  to 
results  from  the  infrared  cameras  (Zappa)  and  thereby  extended  to  smaller  wave  scales.  Breaking 
probabilities,  momentum  flux  and  energy  dissipation  can  be  extracted  from  these  distributions  (see 
Gemmrich  et  al,  2008). 
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Figure  2:  Preliminary  results  on  breaking  crest  length  distributions  A (c)  for  the  Santa  Barbara 
channel  experiment  (left)  and  the  Hawaii  experiment  (right).  Each  curve  is  based  on  a  20  minute 
video  record.  Colour  coding  represents  the  dominant  steepness  KHs  (blue:  high,  red:  low).  The 

dashed  line  indicates  a  c 6 -dependence. 


IMPACT/APPLICATIONS 

This  effort  will  provide  a  far  more  detailed  characterization  of  the  wind  driven  air-sea  interface, 
including  wave  breaking  (whitecaps  and  microscale  breaking).  This  is  needed  to  provide  more 
complete  parameterizations  of  these  processes,  which  will  improve  the  accuracy  of  ocean  optical 
radiative  transfer  models  and  trans-interfacial  image  reconstruction  techniques. 

RELATED  PROJECTS 

The  present  project  is  related  to  the  ONR  project  WAVE  ENERGY  DISSIPATION  AND  THE 
DISTRIBUTION  OF  BREAKING  CRESTS,  (ended  December  2009),  in  which  Andrew  Jessup 
(APL,  UW)  was  the  principal  investigator  and  I  was  a  Co-PI  (via  subcontract).  In  this  project  we 
looked  at  breaking  crest  length  distributions  and  co-located  subsurface  energy  dissipation 
measurements  in  a  strongly  forced  wave  field  in  a  lake  and  in  Pudget  Sound  (Thomson  et  al,  2009; 
Gemmrich,  2010).  While  the  wave  scales  in  RaDyO  and  the  lake/sound  experiments  are  different, 
common  aspects  of  the  data  analysis  have  been  transferred  to  our  RaDyO  data  sets. 
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